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We investigate the time-evolution of a charge qubit subject to quantum telegraph noise produced 
by a single electronic defect level. We obtain results for the time-evolution of the coherence that 
are strikingly different from the usual case of a harmonic oscillator bath (Gaussian noise). When 
the coupling strength crosses a certain temperature-dependent threshold, we observe coherence 
oscillations in the strong-coupling regime. Moreover, we present the time-evolution of the echo 
signal in a spin-echo experiment. Our analysis relies on a numerical evaluation of the exact solution 
for the density matrix of the qubit. 
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The unavoidable coupling of any quantum system to a 
noisy environment leads to decoherence. Understanding 
decoherence is interesting for fundamental reasons (the 
quantum-classical crossover, the measurement problem 
etc.), and is essential for achieving the long dephasing 
times neccessary for building a quantum computer and 
other applications. The paradigmatic models in this field 
(Caldeira-Leggett and spin-boson model jjl |2j [3l [H |5|) 
usually consider a bath of harmonic oscillators. In that 
case, the bath variable coupling to the quantum system 
displays Gaussian-distributed fiuctuations. This feature 
affords considerable technical simplifications, while these 
models are faithful descriptions of real environments like 
the vacuum electromagnetic field or the harmonic crys- 
tal lattice. In other cases (like electronic Nyquist noise 
in a bulk metal), these models represent very good ap- 
proximations. This is a consequence of the central limit 
theorem, applied to the sum of contributions from many 
independent non-Gaussian noise sources. The approxi- 
mation finally breaks down when one couples strongly 
to a few noise sources. This situation is becoming more 
prevalent nowadays, as one studies the coherent dynam- 
ics of nanostructures. The coherence times of solid state 
qubits are often determined by a few fiuctuators tSiL&iSj- 

This challenge has given rise to a number of theoretical 
studies of qubits subject to fiuctuators producing tele- 
graph noise [iiiainiiiaiiaiiiiisiiiainiiiaiTg 

(and other non-Gaussian baths [2Ql [2TT [22l [23l [24l l25]). 

The most straightforward but realistic fully quantum- 
mechanical model consists of a single level tunnel-coupled 
to an electron reservoir [26j. Grishin, Yurkevich and 
Lerner recently studied the long-time limit of this model 
and derived the dephasing rate for a qubit coupled to such 
a fiuctuator [lOj. They found a striking non-analytic de- 
pendence of the dephasing rate on the coupling strength 
and temperature. In the present paper, we take up the 
same model, which may reasonably be termed "quantum 
telegraph noise", now asking for the full time-dependence. 
We find that in the strong-coupling regime (beyond a 
certain threshold) the monotonous decay of the qubit 's 
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Figure 1: (a) Schematic picture of the bistable fiuc- 
tuator: A localized level tunnel-coupled to an electron 
reservoir. (b) Time-evolution of the visibility \D(t))\ 
for classial telegraph noise (top to bottom: v/^ — 
0.2, 0.6, 1.0, 1.4, 1.8, 2.2, 2.6, 3.0). 

coherence turns into temporal oscillations, with complete 
loss of coherence interspersed between coherence revivals. 
We are able to fully include quantum fiuctuations, by a 
numerical evaluation of the exact solution for the quan- 
tum model, and we discuss the behaviour at low tem- 
peratures. We conclude by showing how to extend these 
calculations to spin-echo experiments, relevant for coher- 
ence control. 

Model- We study a single, spin-polarized impurity 
level (Fig. [ija)), tunnel-coupled to a (non-interacting) 
electron reservoir: 

Hb = e^d^d^^^ (^tfc4<^ + h.c.^ + ^^fcc^^fc- (1) 

k k 

Here d) creates an electron on the impurity level of en- 
ergy and tk is the tunneling amplitude to the reservoir 
level /c, of energy Sk (we fix the reservoir's chemical po- 
tential as /i = 0). Below, we always refer to the tunneling 
rate 7 = 27r \tk^ 8{ek — cq). The fiuctuating impurity 
charge Q = d)d couples to a qubit, and the full Hamilto- 
nian is given by (fi = 1, /cb = 1) 

7^ = l^^. + ^Q^.+Ti^B, (2) 



2 



where ^x,y,z ^ire the qubit Pauh operators, A is the 
qubit level spacing, and v is the qubit-fluctuator coupling 
strength. The coupling considered here leads only to pure 
dephasing and not to energy relaxation in the qubit. This 
is a popular and realistic model when discussing the de- 
cay of quantum information during storage. 

We are interested in the full time dynamics of the re- 
duced density-matrix p(t) of the qubit, after preparing it 
in a superposition state and switching on the interaction 
with the fluctuator. Since the interaction, TYint = fQ^^? 
commutes with the qubit Hamiltonian, only the off- 
diagonal elements pij are affected G {T?!})? acquir- 
ing an additional coherence factor D{t)'. 

Pn{t) = Pn{0)e-'^'D{t)- (3) 

Classical telegraph noise. - We first review the clas- 
sical limit for the bath, where the charge Q{t) is a 
stochastic process of the "telegraph noise" type [27J, 
which flips randomly between and 1 (occuring with 
equal prababilities) at a rate 7. This corresponds pre- 
cisely to the high-temperature limit of the quantum 
model discussed here (see below). For a given realiza- 
tion of Q{t), the Schrodinger-equation yields a super- 
position of the qubit 's eigenstates with a random con- 
tribution to the relative phase, (p{t) = —v dt'Q{t'). 
The noise average yields the coherence, D{t) = (e*"^*^^^). 
If the phase were Gaussian distributed, then the co- 
herence would be determined by the variance of (p: 
(e^^(^)) = e^^^(^)>-H^'(^)) . This is not true for clas- 
sical telegraph noise, where the exact result is found to 
be D{t) = e-i(^-^^)^ [cosh((5t) + {j/25) smh{St)], where 
S = — v'^, and 7"^ is the charge correlation time: 

{SQ{t)SQ{0)) = le-^\'\,withSQ{t) = Q{t)-{Q{t)). The 
"interference contrast" of any observable sensitive to the 
relative phase between the qubit 's levels is reduced by the 
factor \D{t)\, which we will term the visibility. Fig. [TJb) 
shows \D{t)\ for different couplings v. Coherence oscil- 
lations appear when v > 7, as (5 becomes imaginary. 
These are qualitatively different from anything observed 
for Gaussian noise, where D{t) cannot cross zero. The 
long-time decay rate of \D{t) \ is equal to ^ (7— ^^7^ — v'^) 
if V ^ 7 and 7/2 if v > 7. 

General exact solution. - In the full quantum model 
[Eqs. ([2]) and ([T])] the coherence can generally [28] be 
written as an overlap, D{t) = (Xb(^)IXb(^))' 
bath-states |xL(^)) IXb(^)) Produced under the ac- 
tion of the qubit being in state |t) or ||). Then the 
coherence is 

D{t) = (^e'i^B--^Q)t^-^inn+^Q)t'^ ^ (4) 

where we average over the thermal state of the electron 
bath. A variety of methods have been applied to calculate 
averages of the form Eq. Q, e.g. linked-cluster expan- 
sions or nonequilibrium Keldysh path-integral techniques 
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Figure 2: Time-evolution of the visibility \D(t)\ for different 
couplings V, for quantum telegraph noise acting on a qubit 
at low temperatures (T/7 = 0.01). The dashed lines show 
the Gaussian approximation. From top to bottom: v/^ = 
0.2, 0.6, 1.0, 1.4, 1.8, 2.2, 2.6, 3.0 (with £0 = 0). 

[Tm [29] . Here we implement a variant of a formula known 
from full-counting statistics [30l [3TJ |32j [33] , which can be 
evaluated numerically efficiently. Given arbitrary single- 
particle operators A, B and C, and their second quan- 
tized counterparts A = XI /c k' ^\'^k'kCk etc., the trace 
tr[e'^e^e^] over the many-body Hilbert space is equal to 
det[l + e^e^e^]. Applying this to Eq. (|4|, we obtain 

D{t) = det (1 - n + e^(^--f ^)^e-^(^- + tO)t^^ 

Here Hb and Q are the single-particle operators cor- 
responding to TLb and Q, and n = /{He) is the 
single-particle equilibrium density matrix, where f{s) = 
{ex.-p{(3e) + 1)~^ is the Fermi-Dirac distribution. This 
formula takes into account exactly the effects of quan- 
tum fluctuations (on top of thermal ones), and the non- 
Markovian features in the fluctuator dynamics that de- 
velop for decreasing temperatures. 

Numerical evaluation.- Our results for the time- 
evolution of the visibility have been obtained by direct 
numerically exact evaluation of Eq. ([5|. To this end, we 
employ a discretization with N equally spaced energy 
levels e e [—W^ W] in a band W These represent 

the single-particle energy eigenlevels of Hb, for which the 
matrix-elements of Q are equal to 

Qap = —JlmG^iu; = Sa)lmG^{u^ = s^). (6) 

Here G^{uj) = {u — Sq -\- ij/2)~^ is the impurity level's 
retarded Green's function and u = N/{2W) is the level 
density. The coherence is obtained by calculating the de- 
terminant of the resulting N x A^-matrix, Eq. Good 
convergence is obtained already for N on the order of 400 
and W = 20. 

Results for the visibility.- In Fig. [2] we show the visibil- 
ity for different couplings v. For small coupling, v/j <C 1, 
the Gaussian approximation works well. It can be ob- 
tained from Eq. by writing det(A) = exp(tr(ln(74))). 
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Figure 3: (color online) Density plot of e^*^*D(t), which is 
real-valued, as a function of time (horizontal) and coupling 
(vertical), (a)-(d) T/7 = 1.0, 0.3, 0.2, 0.01 (with £0 = 0). The 
bold green lines indicate the contours of vanishing coherence, 
D{t) = 0. The dashed line indicates the critical coupling 
strength vl/^i = 1.1, 1.9, 2.2, 2.7 (from (a)-(d)). 

and keeping only the terms up to order v'^ in the expo- 
nent (see also [34J). Equivalently, one may use D{t) that 
would be obtained for a harmonic oscillator bath whose 
two-point correlator is fixed to be (^SQ{t)SQ{0)'^. This 
approximation yields a long-time exponential decay at a 
rate F^^ = v'^ / for T j (agreeing with the results 
for classical telegraph noise, see above). At T = 0, one 
obtains a power law-decay D{t) ^ with an exponent 
a = (4/7r^)('u/7)^, arising from the orthogonality catas- 
trophe [35l[36]. For larger coupling strengths, v/j > 1, 
the Gaussian approximation fails even qualitatively, in- 
dicating the non-Gaussian nature of quantum telegraph 
noise. 

The important feature is the occurrence of visibility 
oscillations beyond a critical coupling strength Vc- The 
visibility vanishes at certain times and shows coherence 
revivals in-between. These features continue to exist in 
the full quantum model. For T ^ 7, it agrees with 
the classical result, where the threshold is = 7. In 
the quantum case (Fig. [2|, we observe a transition to a 
non-monotonous behaviour as a precursor to the visibil- 
ity oscillations, in contrast to the classical limit discussed 
above. Moreover, zeroes in the visibility develop only at 
a larger coupling strength v^, which depends on temper- 
ature T. Another notable feature is the non-monotonous 
evolution of peak heights for v/j > 2.7, unlike the clas- 
sical case. 

To illustrate these points, we have plotted the time- 
evolution of D{t) [excluding a trivial phase factor] as a 
function of the coupling-strength v for various tempera- 
tures [Fig. [3]. At high temperatures, visibility oscilla- 
tions set in at v'^/j ~ 1, whereas for T ^ [Fig. |3ja)] 
the first zero-crossing appears only at v^/7 2.7. 
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Figure 4: Critical coupling strength v^(T) as a function of 
temperature (with £0 = 0,7 = 1). The strong-coupling 
regime is located above the dashed line. At high tempera- 
tures, one has v^{T) 1, according to the classical limit. 



Temperature- dependence of strong- coupling threshold - 
As explained above, the visibility oscillations are a gen- 
uinely non-Gaussian effect. We characterize the on- 
set of the strong coupling regime by the temperature- 
dependent critical coupling v^{T), beyond which the ze- 
roes in D{t) appear. At a fixed temperature T, the crit- 
ical coupling-strength and the corresponding zero in 
D{t) at time are found numerically by a bisection al- 
gorithm. The result is a "phase-diagram" showing the 
critical coupling as a function of T, [Fig. [i]. The 
curve v^{T) separates the v — T-plane into two regions. 
At high temperatures T the critical coupling converges 
to its classical value, ^ 7 [a slight offset in the plot 
is due to limited numerical accuracy]. For low T, it in- 
creases and saturates at a finite value, as D{t;v^T) is 
continuous in the limit T ^ 0, and D{t;v^T = 0) still 
displays oscillations beyond some threshold. This means 
the equilibrium quantum Nyquist noise of the fiuctuator 
is enough to observe visibility oscillations, in contrast to 
the strong-coupling regime studied in [37J, where only 
the nonequilibrium shot noise of discrete electrons could 
yield these effects. 

Spin-Echo - Finally, we investigate the time-evolution 
of the density matrix of the charge qubit in a spin- 
echo experiment, commonly employed to filter out low- 
frequency fiuctuations, whose effect is cancelled in such 
a procedure. Echo protocols were first invented in NMR, 
but they are by now standard in qubit experiments, par- 
ticularly in the solid-state, where they are used to fight 
1/f noise ([38j). At the initial time t' = 0, the qubit 
is prepared in a superposition of its two eigenstates, 
IV^(to)) = 1/72(11) + ID). Then we let the qubit evolve 
according to Eq. ([2]) up to a time t' = t/2, at which we 
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Figure 5: Time evolution of the spin-echo signal, |DEcho(t)| 
(solid line), after applying a vr-pulse at = t/2, in comparison 
with the visibility \D{t)\ for free evolution (dashed line), (a)- 
(d): T/7 = 1.0, 0.5, 0.25, 0.01. {v/'j = 3.0, So = 0). 



perform a 7r-pulse e*^^^ on the qubit, before evolving up 
to time t. Defining U± = ex.p[—i{7iB ± vQ/2)t/2], we 
find the qubit 's final density matrix to be (in analogy 

to Eq. (3)) I^Echo(^) = (ul_UlU-U^y As before, we 
^e this as a determinant in the single-particle 



can rewrr 



Hilbert space, 

^Echo(^) 



det 



(l-n + [/l^|^_t/+n), (7) 



where U± is the single-particle evolution operator. In 
Fig. [5] we compare the echo signal with the free evolu- 
tion. At low temperatures, the fiuctuations are purely 
quantum in origin, yielding a relatively lower weight for 
small frequencies and thus a decrease in the effectiveness 
of the spin echo procedure. 

Conclusion - In conclusion, we have studied the deco- 
herence of a qubit subject to quantum telegraph noise. 
We have calculated the time-evolution of the coherence 
and found a strong-coupling regime with an oscillatory 
time-dependence of the coherence that cannot be mim- 
icked by any Gaussian noise source. We have charac- 
terized this regime via the appearance of the first zero 
in the time-evolution of the coherence and summarized 
the result in a "phase-diagram". Moreover, we have pre- 
sented the time-evolution of the echo-signal in a spin-echo 
experiment and compared it to the coherence. Straight- 
forward extensions of the formulae presented here may 
be applied to discuss the effects of more sophisticated 
pulse sequences [39l HOl HU |l2] which are relevant for 
protecting quantum information storage. 
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